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Abstract 

This paper analyzes the error estimates of the hybridizable discontinuous Galerkin (HDG) 
method for the Helmholtz equation with high wave number in two and three dimensions. The 
approximation piecewise polynomial spaces we deal with are of order p > 1 . Through choosing a 
specific parameter and using the duality argument, it is proved that the HDG method is stable 
without any mesh constraint for any wave number k. By exploiting the stability estimates, the 
dependence of convergence of the HDG method on k, h and p is obtained. Numerical experiments 
are given to verify the theoretical results. 

Key words. Hybridizable discontinuous Galerkin method, Helmholtz equation, high wave 
number, error estimates 

1 Introduction 

The numerical solutions of Helmholtz problems have been an area of active research for almost 
half of a century. Because of the well known pollution effect, the standard Galerkin finite element 
methods can maintain a desired level of accuracy only if the mesh resolution is also appropriately 
increased. In order to remedy this problem and to obtain more stable and accurate approximation, 
numerous nonstandard methods have been proposed recently (cf. [21] ). One type of methods applies 
the stabilized discrete variational form to approximate the Helmholtz equation, which includes 
Galerkin- least-squares finite element methods [10.19.25j, quasi-stabilized finite element methods [5], 
absolutely stable discontinuous Galerkin (DG) methods |14H16| and continuous interior penalty 
finite element methods (CIP-FEM) [29|. Other approaches include the partition of unity finite 
element methods [5j,l24|l2B]. the ultra weak variational formulation [9|, plane wave DG methods 
[2)120] . spectral methods [27], generalized Galer kin/finite element methods [71123]. meshless methods 
[6], and the geometrical optics approach [13]. 

Discontinuous Galerkin methods have several attractive features compared with conforming 
finite element methods. For example, the polynomial degrees can be different from element to 
element, and they work well on arbitrary meshes. For the Helmholtz equation, the interior penalty 
discontinuous Galerkin methods (cf. [21 [15]) and the local discontinuous Galerkin methods [16] 
perform much better than the standard finite element methods, and they are well posed without 
any mesh constraint. Despite all these advantages, the dimension of the approximation DG space 
is much larger than the dimension of the corresponding classical conforming space. 
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The hybridizable discontinuous Galerkin methods were recently introduced to try to address 
this issue. The HDG methods retain the advantages of the standard DG methods and result in a 
significant reduced degrees of freedom. New variables on the boundary of elements are introduced 
such that the solution inside each element can be computed in terms of them. In particular, 
element by element, volume degrees of freedom can be parameterized by the surface degrees, and 
the resulting algebraic system is only due to the unknowns on the skeleton of the mesh. For a 
comprehensive understanding of the HDG methods, we can refer to [IT] for a unified framework 
for second order elliptic problems and to [22J for the implementations. 

In |17| . the authors give error estimates of the HDG method for the interior Dirichlet problem 
for the Helmholtz equation, but it is under the condition that C K «;/iM™ m and C K K/iM™ ax are 
sufficiently small, where C K is a constant which is dependent on k but is not characterized explicitly, 
and M™ n , M™ ax depend on the parameters defined in the numerical fluxes. Motivated by this work, 
the primary objective of this paper is to analyze the explicit dependence of convergence of HDG 
method for the Helmholtz equation on k, h and p. In this paper, we consider the Helmholtz equation 
with Robin boundary condition which is the first order approximation of the radiation condition: 

-Au-K 2 u = f in ft, (1.1) 
du 

— h iku = g on dQ, (1-2) 

on 

where O C R d , d = 2, 3 is a polygonal/polyhedral domain, n > is known as the wave number, 
i = v 7 — 1 denotes the imaginary unit, and n denotes the unit outward normal to Oil.. 

The main difficulty of analyzing the Helmholtz equation lies in the strong indefiniteness of the 
problem which makes it hard to establish the stability for the numerical approximation. For the 
HDG method, we use a duality argument to obtain the stability estimates of the numerical solution. 
In the analysis, a crucial step lies in the derivation of the dependence of convergence on p. We 
utilize the explicit error estimates of L 2 projection operator (see Lemma 14. 4p to overcome this 
problem. Then we obtain that the HDG method for the Helmholtz problem (|l.ip - (|1.2p attains a 
unique solution for any k > 0, h > 0. Furthermore, the stability results not only guarantee the 
well-posedness of the HDG method but also play a key role in the derivation of the error estimates. 

The duality argument can not be directly applied to establish the error estimates. Thus, we first 
construct an auxiliary problem and show its HDG error estimates by the duality technique. Then, 
combining the stability estimates, the error estimates of HDG scheme for the original Helmholtz 
problem (jl.ip - (|1.2p are deduced. Let Uh and q h be the HDG approximation to u and q := iVu/n 
respectively. We obtain the following results: 

(i) The following stability and error estimates hold without any constraint: 

(K/^ lip 1 \ / Kj 2 ]% \ 

i + —^r) ll/llo,n + [ 1 + —^-) llflkan, 

. /nh 2 K 2 h 2 k 5 /i 4 \ , r . ~ . 

\\U - U h \\ ,Q < —t- H 2~ H 4~ ) M if'9), 

\ £r p A / 

. (Kb. ^h 2 k 6 /i 4 \ ^ „, 
K||«-9 fc lkn< [— + -^ r + -^-)M{f,g), 

where / := —if/n, g := —ig/n and M(f,g) := ||/||o,n + ||<?||o,dfJ- We use notations A < B and 
A > B for the inequalities A < CB and A > CB, where C is a positive number independent of 
the mesh size, polynomial degree and wave number k, but the value of which can take on different 
values in different occurrences. 
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3 1,2 

(ii) Suppose % 1) there hold the following improved results: 

IWIo.fi + ||9hllo,fi % HZ llo.fi + Ibllo.Sfii 
„ , (nh 2 n 2 h 2 \ ^ r ,~ „. 

\\U - UftHo.fi < [—J- + ~ §~ ) M (f'9), 

— + — T -)M(f,g). 
p p A / 

Comparing to the estimates in /ip-IPDG method for Helmholtz problem, we find that the condition 

3,2 

for the above improved results weakens the mesh condition < 1 which is requested in |15| . For 

3 ul 

the estimates under the mesh condition ^-p- > 1, the results in (i) can not be directly applied, but 
we may still get the following improved estimates. 

(iii) Suppose ^j- > 1, there hold 

K 2 h 2 

\\u - UftHo.fi < — =-Af(/, 5), 

«l|flf " 9fcl|o,n < (— + ^r)M(f,g). 
p p A ' 

We remark that in this work the local stabilization parameter to determine the numerical flux 
in the HDG scheme is always selected as r = -0- (see (|6.8p ). Our numerical results show that the 
predicted convergence rates are observed. 

The organization of the paper is as follows: We precisely define the HDG method for the 
Helmholtz equation and give some notations in the next section. Section 3 is dedicated to the 
characterization of the surface degrees Uft. In section 4, we derive the stability estimates of the 
HDG method. The error estimates of the auxiliary problem are carried out in section 5 while 
section 6 states the main results of this paper, i.e., the error estimates of the HDG method for the 
Helmholtz equation. In the final section, we give some numerical results to confirm our theoretical 
analysis. 

2 The hybridizable discontinuous Galerkin method 

The HDG scheme is based on a first order formulation of the above Helmholtz equation (jl.ip - (jl.2p 
which can be rewritten in mixed form as finding (q, u) such that 

iKq + Vu = in ft, (2.1) 
iKU + divq = f in ft, (2-2) 
— q • n + u = g on 5ft. (2-3) 



Existence and uniqueness of solutions to ()2.1|) - ()2.3p is well known and it is proved in [12] that they 
satisfy the following regularity result: 

« _1 ||u||2,fi + ||u||i,n + lklli.fi < M(f,g). (2.4) 

We consider a subdivision of ft into a finite element mesh of shape regular triangle T in M? (or 
tetrahedron T in M 3 ) and denote the collection of triangles (tetrahedra) by 7ft, the collection of 
edges (faces) by £ft, while the collection of interior edges (faces) by £® and the collection of element 
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boundaries by dTh ■= {dT\T E Th}- Throughout this paper we use the standard notations and 
definitions for Sobolev spaces (see, e.g., Adams PQ). 

On each element T and each edge (face) F, we define the local spaces of polynomials of degree 
p > 1: 

V(T) := (V p (T)) d , W(T) := V P (T), M(F) := V P (F), 

where V P (S) denotes the space of polynomials of total degree at most p on S. The corresponding 
global finite element spaces are given by 

-2 

-2 



V p : = {v E L 2 (tt) | v\ T € V(T) for all T E %}, 



< 



{w E I/(Q) I w\ T E W(T) for all T E %}, 
{H E L 2 (£ h ) | (jl\f E M(F) for all FeS h }, 

where Z 2 (ft) := {L 2 (fl)) d and L 2 (f h ) := U F&£h L 2 (F). On these spaces we define the bilinear forms 

{v,w) Th := (v,w) T , {v,w) Th ■= ( v ' w )t, an d (v,w) 9Th := ^ (v,w) dT , 
Te% TeT h TeT h 

with (v,w)t '■= j T v -wdx, (v,w)t '■= J T vwdx and (v,w)qt '■= j aT vwds. 

The hybridizable discontinuous Galerkin method yields finite element approximations (q h , Uh,Uh) E 
V p h xW%x Ml which satisfy 



{™q h ,r) Th - (u h ,divr) Th + (u h ,r ■ n) dTh = 0, (2.5) 

(iKu h ,w) Th ~ {Qh^ w )r h + (Qh ■ n ,w)dT h = (f,w)r h , (2.6) 
{-Qh • n + u h ,-p) d n = (g,~p)an, (2.7) 
{q h -n,]l) drh \ dn = Q, (2.8) 

for all r E V p h , w E W^, and fi E M^, where the overbar denotes complex conjugation. The 
numerical flux q h is given by 

Qh = Qh + T ( u h ~ u h )n on dT h , (2.9) 

where the parameter r is the so-called local stabilization parameter which has an important effect on 
both the stability of the solution and the accuracy of the HDG scheme. We always choose t = Sc 
in this paper. The error analysis is based on projection operators which are defined as follows 

n h : L 2 (n) ^V p h U h : L 2 (Q) w p 

for any T E Th, they satisfy 

(U h q, v) T = (q, v) T for all v E V(T), (2.10) 

(U h u,w) T = (u,w) T for all w E W(T). (2.11) 

We conclude the introduction by setting some notations used throughout this paper. Let the broken 
space H l {£lh) be defined by 

H\Sl h ) := {v : v\ T E H\T), VT E T h }, 

the seminorm of which is 



\ v \l,a h 



TeT h 
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The trace of functions in H 1 ^) belong to T(T) := U TeTh L 2 (dT). For any <f> G T(T), and 
v G (T(r)) d , if e € £JJ, e = 5T+ n ST - , we set 

{{0}} := ^(</> + + <T), M := + • n+ + <j>~ ■ n 

and 

{{i>}} := ^(f" 1 " + v ~)i l v i := v+ ' n+ + u_ ' n • 
For e G 97h n dU, we define 

M := 0, M := • n, {{v}} := «, [«] := v ■ n. 

3 The characterization of Uh 

One of the advantages of hybridizable discontinuous Galerkin methods is the elimination of both 
q h and Uh from the equation and obtain a formulation in terms of ■% only. In this section, we show 
that Uh can be characterized by a simple weak formulation in which none of the other variables 
appear. 

First we define the discrete solutions of the local problems, for each function A G M?, {QxMx) £ 
V(T) x W(T) satisfies the following formulation 

(iKQ x ,r) T -{U x ,di^)T = -(\r T n)dT for all r G V(T), (3.1) 
(iKW A ,^)T-(QA,V^)r + (QA-n,^)oT = for all w G W(T), (3.2) 

where Qx ■ n = Qx ■ n + r(W A - A). For / G L 2 (Q), (Qf,Uf) G V"(T) x is defined as follows 

(iKQ f ,r) T - (U f ,divr)T = for all r G V(T), (3.3) 
(i^f,w) T -(Qf,Vw) T + (Qf-n,w) aT = (f,w)T for all w G W(T), (3.4) 

where Qf ■ n = Qf • n + rUf. Next we show that the local problem (|3.1|) - (|3.2|) is well posed. The 
uniqueness of (|3.3[) - (j3.4[) can be deduced similarly. 

Lemma 3.1. There exist a unique solution {QxMx) G V(T) x VF(T) to the local problem ()3.1j) - f)3.2[) . 

Proof. Since it is a square system, to prove the existence and uniqueness of its solution, it is enough 
to show that if A = 0, we have that Q\ = 0, U\ = 0. Taking r = Qx in f)3. 1 [) . w = U\ in (|3.2p . we 

get 

\k{U x Mx)t - 'ik(Q x , Qx)t + t{U x Mx) 9 t = 0, 

which means U\\qt = 0. Back to (|3.ip . we derive that 

-ikQx = VU X - 

Inserting the above expression into (|3.2p . 

AU X = -K 2 U X 

is deduced, which implies that Ux = 0, Qx = 0. □ 



5 



It is worth noting that the solution (q h ,Uh) in (|2.5|) - (|2.8|) is exactly correspond to the following 
relationship 

q h = Qu h + Qf, u h = Uu h +u f . 

And Uh is the solution of the following formulation 

ah(u h ,fi) = bhip) for all \i G M%, 

where 

a h (X,fi) : = -(lQ\},~p)dT h + 
h(l*) : = (lQfhJI)dT h + (0,P)an- 

4 The stability of the hybridizable discontinuous Galerkin method 

The goal of this section is to derive stability estimates. We first cite the following lemma which 
provides some approximation results that will play an important role later. A proof of the lemma 
can be found in [Hl2"8"]. 

Lemma 4.1. Let T be a standard square or triangle. Then there exists an operator tt p : i/ 1 (T) — > 
V P (T) such that for any u G H l (T) 

\\u-TT p u\\^f <P _1 |u|l,T- 

Moreover, if u G H 2 (f), 

\\u-TT p u\\ 0jf <p~ 2 \u\ 2jf , (4.2) 

\u-7T p u\ hf <V' l \u\ 2 ,f- ( 4 - 3 ) 

Using the standard scaling technique, we can get the following approximation results. 

Lemma 4.2. For any T G 7~h, there exists an operator 7r^ : H l (T) — > V P (T) such that for any 
u G iJ x (r) there holds 

ll u - ^llo.T < -|«|i,t, (4.4) 
p 

Moreover, if u G H 2 (T), 

\\ u - ^h u h,T ^ (^) kb.T, (4.5) 

| ii - 7T?u|l,T < -|li| 2 ,T- (4.6) 

p 

We also need the following trace inequality and refer to [28] for the proof. 
Lemma 4.3. For any T G Th and i> G V P {T) 

[Mlo,9T < p/i~ 5 lb||o,T- 

Now we derive the following approximation properties of the projection operator 11^ which is 
defined in (|2.1ip . For the sake of simplicity, the proof is restricted to 2-d case. 
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Lemma 4.4. For any T G l~h, u G H l (T), the projection operator II/j satisfies 



| it - n^nllo; 



< _ 

rsj ! 

P 



\u - H.hV\\a,dT 



< 



Moreover, if u G i? 2 (T) 



\u - n h it||o,3 



< 



h\2 



«1,T- 



P 



F 2,T, 



| it - n h u\\ ,dT < y~) \(i\2.:r- 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



Proof. An important property of L 2 (VL) projection operator II/j is that 

\\u - U h u\\o :T < mi veW (r)\\u ~ v\\o,t, 

hence (|4.4p and (|4.5f) imply (j4.7j) and (|4.9p . Let Ti and T 2 be the standard triangles with linear 
mappings x = F±(x) and x = ^(x) respectively, see Figure Q] for illustration. For any v G W(T), 
define v(x) := v o i*i(x) and -u(x) :=»o ^(x). 




Figure 1: The triangle T and its reference triangles T\ and T%. 
For any -u(x) G i/ 1 (Ti), we have 



|£>(x)| = v(x) ■ v(x) = v(xi, 0) • v(xi,0) + 



■t'2 



d(v(xi,rj) ■ v(xi,rj)) 



drj 



\v(xi,0)\ +2Re / u(xi,7?) 



X' 2 



9?5(xi, rj) 



drj 



drj. 



Hence 



r%2 

|i)(xi,0)| 2 = |-0(xi,x 2 )| 2 - 2Re / v(xi,rj) 

Jo 



dv{x\,rj) 



drj 



drj. 
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By using integration on T\ we can deduce 

1 pi rl—Xi 

\v(x\, 0) | 2 ( 1 — xi)dx\ = I I \v(xi, 0)\ 2 dx2<ixi 



o Jo 



nl-xi f-x 2 dv(xi,i]) 
/ v(xi,i]) — - — dr]dx2(lxi 
Jo 



\v{x\, X2)\ 2 dx2dx\ — 2Re / / / v(xi,rj) — - — d^drjdxi 



Ti JO JO 



1 >' x Xl dv(xi,rj) 



/ \v(xi,X2)\ 2 dx2dxi - 2Re / / (1 — x\ — r])v(xi, rj) 
Jfi Jo Jo 



di] 



di]dx\ 



dv (xi , X>2 ) 

\v(xi, X2)\ l dx2dx\ — 2Re / (1 — x\ — x<i)v{x\, X2) ttt dx^dxi. 

0x2 



T, 



Take v = u — ILu, where Ilrt G V p (Ti) satisfies (ILu, w)f = (u,w)f , \/w G V P (T\), and note that 

ITI ■ — ■ 

(u - IL h U, V) T = pr- (u - IlftU, 0)^ = 0, 
l T l| 

which means ftu = D^u. By Lemma 14.21 and scaling technique we get 
((u — ILu)(xi, 0)) 2 (1 — x\)dx\ 



and 



f f d(u — ITn) 

= / (u — Tlu) 2 dx2dxi — 2Re / (u — ILu) (1 — x\ — X2)dx2dx\ 

Jfj Jn 0x2 

= [ (u — tiu) 2 dx2dxi — 2Re f (u — ILu) - ~„ p "^ (1 — &\ — X2)dx2dx\ (4-H) 

Jt x Jt x (JX2 

< ||u-n«[|jj^ + 2\\u - iiu\\^ fi \u - Tr p u\ ltfi 

h 2 

< h~ 2 \\u - U h u\\o,T + h \\u - Yl h u\\o tT \u - n p h u\x,T < ^3 Mi,T> 



((u — LLu)(x\, 0)) 2 (1 — x\)dx\ 



J (u — Tlu) 2 dx2dxi — 2Re J (u — — x\ — X2)dx2dx\ ^ 

|2 



< ||u - nu|| 0)fi + ||n-nu|| 0jfi |-u| ljfi 

< hT 2 \\u - LL h u\\l T + h~ l \\u - n h it||o,T|«|i,T < v~ 1 W\\t- 



Now we map T to T2 and similarly we can derive 

,2 



and 



\ ((u -nu^x^orfxidx! < —\u\l T , (4.13) 
J ((u - ILu)(x 1, 0)) 2 Xl dx x < p-Min (4-14) 



where flu £ V P (T2) satisfies (TLu,w)f = (u,w)f , Vtt) G V P (T2) and ITu = ILj-u. Since v{x) = v(x) 
on e, summing up (|4.1ip . (|4.13[) . and (|4.12p . (|4.14p respectively and noting that e is not particularly 
chosen, the lemma is proved. □ 

Remark 4.1. In this paper, we only deal with meshes consisting of triangles or tetrahedra, but 
we should note that Lemma 14.41 can be extended to the meshes constituted with rectangles or 
hexahedra. The proof is similar with Lemma 14.41 and can also be found in [8] . 

Lemma 4.5. Let (q h ,u h ,u h ) £ V p h x W% x M v h be the solutions of (&J^-$ZM). There hold 

i~\\u h - u h \\l,aT h ^ ll/MIKHo,fi + ||p|lo,9fi> ( 4 -15) 
K hh\\o,n ^ K IMIo,n + ll/llo,n + IMItW ( 4 -16) 

Proof. We choose r = q h , w = Uh, ^ = Uh in (|2.5p - (|2.8p and get 



(™q h ,q h )T h ~ (u h ,divq h ) Th + (u h ,q h -n) dTh = 0, (4.17) 

{\KUh + div q h ,Uh)% + ( T (uh - u h ),Uh)ar h = [f,u h ) Th , (4.18) 

(Qh • n > u h)dT h = ( u h, u h ) dn - (g, u h ) 9n . (4.19) 

Using (|4.19p . the complex conjugation of (|4.17p can be rewritten as 

-( iK Qh:Qh)T h ~ (div<z h ,u^)r fe - (r(u h - u h ),u h ) aTh + (uh,Uh)dn = (g,Uh)an- (4.20) 
Adding (|4.2U[) and (|4.18|) together, the following identity is deduced 



(mu h , u h ) Th - (i/eq/i, q h )r h + (r(u h - u h ), (u h - u h )) aTh + (%, u h )an 

= (f,uh)r h + {g,uh)an, 

which implies the lemma. □ 

Next we use a duality argument to estimate the stability of Uh- Given Uh £ L 2 (Q), we introduce 
the dual problem 

-m* + W = inO, (4.21) 
div*-iK* = u /l in CI, (4.22) 
* • n = * onffi. (4.23) 

In the following lemma, we give some explicit bounds for ^ and 3>. 
Lemma 4.6. For and defined above, they admit the following estimate: 

ll*Ho,n + « -2 ||*lkn + Ol^lkfi + II* llo,en + < IKM- (4.24) 

Proof. In fact, ^ satisfies the following equation 

+ k 2 ^ = \Ku h in VL 

• n = ik^ on dfl. 

In [12], it is proved that 

\\^\\o,n + K~ 2 ||*|| 2 ,n + o < ||ith||o,n- 
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Since \P satisfies the following weak formulation a(^f,v) = (iKUh,v), where 

a(^,v) : = -(V*,W) + k 2 (^,v) + iK(y,v) dn . 
Testing the above formulation by v = ^ and taking the imaginary part yields 

K ll*llo,9Q < «l|«fcl|o,n||*||o,n < «ll«hllo,n> 

which finishes the proof of this lemma. □ 

Now we are ready to derive the stability of Uh, which plays an important role in the error 
analysis for the Helmholtz equation. 

Theorem 4.1. Let (q h ,u h ,u h ) G V P h x W% x be the solutions of (pT5|) - (pI%]) . Then 

Kllo,n < (l + ^V) ||/||o,n + (l + h\\ o,dn, (4-25) 

\\q h \\o,n < (l + [|/[|o,n + (l + ^y) \\gh,d€i, (4.26) 

< ({—y+ph-i) ((1 + -^-)||/|| 0l n + (1 + — )Nlo,an). (4.27) 
Proof. Using (|4.22p . Green formulation and the definition of projection operators, we obtain 

(u h ,u h ) Th = (u h ,div& - iK^!) Th = (u h ,div®) Th +in(u h , ^) Th 
= -(Vu h ,l>)T h + {u h ,W- n) 9Th + i/c(« h ,*) 

= -(V%,n h *) r , + (n h ,n h $ • n) aTh + {u h , ($ - • n) dTh + iK(«/„ W)tj; 

= (u h ,divU h ^) Th + (u ft , (* - Il h 3>) ■ n) dTh + i/c(« ft ,II h *)7^. 

Hence using (|2.5p and the fact that $ • n is continuous across the inner edges, the above equality 
can be rewritten as 

(u h ,u h ) Th = (iKq h ,Il h $>) Th + (u h ,IL h & • n) dTh - (u h ,$ ■ n) dTh + (u h , $ • n} 9Th 
+ (* - n h *) • n) aTh + i/c(it h ,n fc *) Th 

= (™q h ,&)% + - (* - n />*) • n )sr h + • «}an + wK,IIft% r 
Green formulation and (|4.2ip indicate that 

(iKq h ,*)r h = (q h ,-V^)r h = (divq^F)^ - (qr h • n,"^) aT;i = (divqf A ,n^*)r h ~ (Qh ' »»*W 

Combining (gSD-fgS]) and fl03|) gives 

(uh,u h )r h = {iKu h + div q h ,U h ^) Th - (q h ■ n,^) dTh 
+ {u h - u h , ($ - n h $) • n) dTh + * • n) m 

= {f,TTi^)r h - (r(u h - u h ),Th^) dTh - (q h ■ n,y) 9Th + (q h ■ n,^) dTh 

- (q h ■ n, ty) dn + (u h - u h , (* - 11^*) • n) dTh + (u h , $ • n) dn 

= {f,TTi^) Th + {r(u h - -n h V) dTh + (u h - u h , ($ - n h $) • n)ar h + (P M #, 
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where we have used that the normal component of q h across interelement boundaries is continuous 
and P M g G V v (dT h n <90), (Pm9, v)an = {9, v)au, G V p {d% n dSl). So we can get 

Kllo.fi ^ ||/||o,n||*||o,n + r||(tt A - u/OlkarJ* - W\\o,dT h 
+ \Wh - Uh\\o,dT h \\® ~ flh®\\o,ar h + l|5l|o,9fi||*||o,5fi- 
Applying Lemmas I4.4H4.5I and the regularity estimate (|4.24jl , we get 

|2 < 



|p/illo,fi ~ ll/l|o,fi|p?i||o,fi + llfl'llo.anllw/illo.n 

+ ( r ^ 2 (~) f +' r "^(^)')(ll/llo I nll"ftllo,n + lbllo,sn)Kllo,n. 

Note that we choose r = to get the minimum of the term t5k 2 (^ + t _ 5k(^)2. Eliminating 
|| n /i||o,fi from both sides of the equation, we can get 

3 h I i 

IKIIo.fi < C||/||o,fi + CIMIo.afi + ^ 5 -(ll/llo,filKllo,fi + llfllo.an) 

< C(l + ^)||/|| ,fi + C(l + — )||5||o,8fi + SIKIIo.fi- 

Choosing 6 < |, (|4.25p is obtained. Using (|4.16p . the bound for q h is deduced. According to 
Lemma 14.31 

IWIo.dfi < ph~ 2 1 1 u ft || 0) n, 
which combined with (|4.25j) . (|4.15p and the triangle inequality yields (|4.27p . □ 

5 Error estimates of an auxiliary problem 

In this section, we derive the error estimates of the solutions for the auxiliary problem 

ikQ + W = in n, 

div Q — ikU = f — 2\ku in 
— Q ■ n + U = g on <9$7, 

where u, f and g are determined by the problem (|2.ip - (|2,3p . The HDG scheme of this problem is 
to find (Q h , C/ft, Aft) G V p h x W% x such that 



(inQ h ,r) Th - (U h ,dWr) Th + (X h ,l^n) aTh = 0, (5.1) 
-{\KU h ,w) Th - (Q h ,Vw) Th + (Q h • n,w) dTh = (f - 2iKu,w) Th , (5.2) 

(-Qh ■ n + K,~P)m = {g,~P)an, (5.3) 
{Q h -n,-p) dTh \ d n = 0, (5.4) 

for all r £ V?, u; G and /i G Mf, where 

Qh = Qh + T ( u h - Ah)ri on dTh- (5.5) 
Inserting the expression of Q h into Q5.3P and (|5.4p . we obtain that on the edge e G dTh\d£l 

^ = {{U h }} + ^lQhl 
Q h -n = {{Q h }}-n+ T -{U h }n, 
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and on the boundary edge e E dTh H 90 



Q h -n + rU h P M g 
1 + T 1 + r 

1 + T 

We can substitute the above expressions into (|5,ip - ([5.4p . and get the equivalent formulations of Q h 
and Uh as follows: 

B 1 (Q h , U h ;r,w) := {^Q h ,r) Th - (U h ,db^) Th + £ {({{U h }}, [r]) e + {^[QJ, [r]) e ) 



2r' 

h 



+ (-r- ,r • n) dn + (— — f/^,?^ n) dn = -(— — #,r • 

1 + T 1 + T 1 + T 

B 2 (Q h ,U h ;r,w) := -(iKU h ,w) Th + (divQ h ,w) Th - £ (IQJ,W) B + £ {^M; N 



(5.6) 



(5.7) 



"Van' 



for allr £V p h ,w £W%. Define 



A(Q h ,U h ;r,w) := Bi(Q h ,U h ;r,w) + B 2 (Q h ,U h ;r,w), 

and 

^iCQa)^;^™) : = Bi(Qh, u h\r,w) + B 2 (Q h ,U h ;r,w). 
An obvious observation is that 

>l(O ft ,^;g A ,^) = i«(||g A ||g in + ||i7 ft ||g >n ) + X; (^rlllQJIIo,e + JlMllU 



1 r 

+ r-r-ll^llo,«i- 

1 + r 1 + T 



(5.8) 



Since the formulation is consistent, we have 

A{<1 -Qhi u ~ u h.; r, w) =Ai(q- Q h , u - U h ; r, w) = 0, (5.9) 

for all r E V?, E W|\ We first give the error estimation of the flux Q h , and then use the duality 
argument to bound the L 2 -error of the discrete solution Uh- 

Theorem 5.1. Let Q h and Uh be the solution of (|5.6p - (|5.7p . Denote e q := q — Qh = Q— H^q+II^eq, 

e u := u - Uh = u - U h u + II fe e n and 

e := «||n h e,||g )n + K||n ftett ||g )n + ^(^||[n ft ej|^ e + ^||pl h e u ]||g ie ) 



+ T—\\ u he g • n\\ om + — — ||n h e u || m - 

1+T 1+T 



Then there hold the following estimates: 



E<^-M 2 (f,g), (5.10) 

\Q-Qh\\o,a<-M(f,g). (5.11) 
P 
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Proof. Direct calculation shows that 

Y / Z®-nds=J2 [UU®}}ds+J2 fUM®}ds, (5.12) 



TeT h JdT ee£ h Je ee s2 



h 



for all G (H 1 (n h )) d and £ G H 1 ^). According to (JESJ) and (i5J2]l we have 

^(n fe e q , n/je u ; n fc e g , Il h e u ) = A(U h q - Q h ,U h u - U h ; U h e q , U h e u ) 
= A(U h q - q, U h u - u; U h e q , U h e u ) + A(q -Q h ,u- U h ; U h e q , U h e u ) 
= A(U h q - q, Il h u - u; Il h e q , U h e u ) 

= Y (({{U h u - u}}, [n^]) e + (^lU h q - qj, [TT^,, , 



i (Tlhq — q) n — . , r — > 

+ ( — ' n * e « " n /dn + ( Y+^^ hU ~ " n /an 

+ m h eul{{U h q-q}}) e + £ <^pl fc e u ], [II^u - u\) e 



+ r— — {n h e u ,Il h u - u) dn + — — (U h e u , (Il h q - q) ■ n) 9Q . 

1 + T 1 + T 

Using (|5.8p and the Young's inequality we obtain 

\E < C s t Y \\Tl h u - u\\l e + 5 Y, ^m h e q \\\l e + C s ± Y \\U h q - qf 0je 

+ -r^-\\U h q - q\\l dn + __||n h e 9 • n\\l m + C S T\\U h u - u\\l m 
1 + r 1 + r 

+ c 5 r- 1 Y \\nhq-q\\ 2 ,e + $Y JlP^JIIo^ + ^^rlln^llo.en 

^ I ITT l|2 

+ 7aT7) l|n ^- 9ll °^- 



Choosing 5 < j an d taking advantage of Lemma \4A\ we get 

\e<tY Pfcu - «||g, e + r- 1 ^ ||n fc g - g||g >e < ^M 2 (/,5), 

and hence ()5.10j) is deduced. Then (|5.1ip follows from Lemma 14.41 and ()5.10p . □ 

Next we establish an error estimate for Uh , we perform an analogue of the Aubin-Nitsche duality 
argument to get the convergence rate. First we begin by introducing the dual problem 

-i/t* + W = in O, 
div + ik^ = e u in $7, 

$ . n = ^ on <9J7, 

and prove its regularity estimations. 

Lemma 5.1. Let and \l/ be defined above, then they admit the following estimate 

||*||i,n + ||*||i,n + || 2 ,n < ||e u || ,n- (5.13) 
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Proof. Direct calculation shows that satisfies the equation as follows 

- k 2 $> = \ne u in 
X7q> ■ n = ik^ on dQ. 

It is well known that ^ is the solution of the following weak problem, for all v G H 1 ($7) 

a(^,v) := (V*, W) + k 2 (^,v) - iK{^,v) dn = (-iKe u ,v). 

Taking v = *S>, we get 

ll^llo.n < K _1 ||e u || ,Q, 

and 

ll*lli,n £ l*li,n + ll*llo,an £ «l|e«||o,n||*||o,n < Ikllo.n, 

where we have used Poincare inequality. The regularity theory for the Laplace problem (see Chap 
2 of [18] ) gives the bound for |^| 2,n, 

\^k,n < ||i«e u - « 2 *||o,n + ||iK*||^i 

< K||e u ||o,n + «||*||i,n ^ «||e«||o,n- 

Combining the definition of 3? and the above estimates completes the proof of (|5.13p . □ 

Now for any 0,r G (i7 1 (Q/i)) rf and G we define the following bilinear form 

A(®,C;r,w) := -(i«e,r) Th - (e,di^F)r h + E «€ej,[r]) e + (^;[e],[r])e) 



+ ( rj — , r • n) fln + < — — e, r • n) 9Q + u;) Th + (div 0, w) Th 



• n . , r 

T+7' r " n ' sn + 'iT'r 

E ([»]. W). + E <^K1. I^>e + (iT7^ - e • n),iu 



Direct calculation shows that 



A(®,Z;r,w) = Ai(-r,w;-®,0- (5.14) 

Moreover, the consistency of the bilinear form implies that 

A(&,^;-e q ,e u ) = (e u ,e^), (5.15) 

where e q ,e u and <&,\P are defined in Theorem 15.11 and Lemma 15. 1 1 respectively. 
Theorem 5.2. Let q fc and £7^ be the solution of (|5.6p - (|5.7p . There holds 

||u-^||o,n<^M(/,fl|). (5.16) 
Proof. Using (joTTB]) . (joTTS]) and ([5U]) . we have 



-«llo,fi 



*; -e q , e u ) = Ai(e q , e u ; *) 



= Ai(e q , e u ; -* + n fc *, * - n h tf ) = ^(I^e,, n fc e u ; + 11^, * - 11^) (5.17) 
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Denote 

T x := Ai{U h e q ,H h e u --^ + U h ^^ -Ii h ^!) 

and 

T 2 := Ai(q - U h q, u - U h u; -* + n h *, - U h ^>). 

Then we estimate the above two terms respectively. By (|5.12|) and the property of the projection 
operators we can rewrite T\ as 

T x = - (PW, {{Th& - $}}) e + E (^-plfce,], [n fc * - *l) e 



2r' 

1 



— — (n h e u , (n h * - *) • n) an + ^— — (n^e, • n, (11^* - *) • n) 
1 + r 1 + t 

r 



/an 

- T 

2 (pi fte ,l, {{^-n^}}) e + E ^([n, e j, |[^-n h *]) t 



1 + r 1 + r 

Applying the Cauchy-Schwarz inequality, we have 



|2i| < ( E T\m h e u }\\l e y ■ (r- 1 E H n fc* " * 



i i 

2 \ 2 

0,e ; 



1,1 .1 

2 



2t ,iil ^„u,ey y 2r 



l 1 
|2 \2 / -lMTT./R *l|2 \ 2 



i i 

r lll"Ae,j|||o, e j a • (r E ll*-M||LV 
i i 



+ ( E 7llPifcc,]||g ie ) 3 • (r E II* - n fc *||§ 
+ (E pI n ^JIIo ie )" • g E ll*-Hfc*llS,«)* 



+ (y^P^IIiUj)* • (t^H* - n^HiW) 5 

Then the upper bound for T\ follows from Lemma 14.44 Theorem 1 5 . 1 1 and the regularity estimation 
([5TT3]) that 

|r 1 |<^M(/,#||e u ||o,n. (5-18) 
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Similarly we use the property of the projection operators and get 
T 2 = iK(q - U h q, Il h 3> - *) Th + (V(n - vr£u), Il h 3> - ®) Th 

- E (I« - n ft «],-pi fc *-*a> e + X) - n h «zj, [n fc * - *j) 



1 



+ z : — ((q - U h q) ■ n, (Il h $ - *) • n) dn + — — ((u - U h u), (Il h $ - *) • n)an 
1 + r 1 + r 



m{u - Ii h u, - n h *) Th - (q - U h q, V(tf - <*)) Th 

r 
2* 



+ X («9 - n^}}, I* - n h *]) e + J2 - n h «], - n^]) e 



+ — — (« - n fc u, * - n h *)an - ^— ((<? - n^q) ■ n, (* - iL^))^. 

1 + r 1 + t 

Hence 

|T 2 | < n\\q - U h q\\ 0i a\\U h ^ - *|| ,n + |u - 7r£u|i A J|n h $ - *|| ,n 



+ ( X ll« - n*«||g, e ) 2 • ( X ||n h * - *||g ie J 2 + 1| 9 - n fc g|| ,n|* - <*|i,n h 

e££ h e££ h 



i 



+ (X>" n^llo.e) 2 • ( E II* " n ^Ho,e) 2 + K W U ~ n ^ll0,n||* " 

e££ h ee£ h 



1 1 

2 



+ r(Ell«-nHlo >e ) 2 -(Ell*- n ^Ho, 

e££ h ee£ h 

Using Lemma H31 Theorem 15.11 and (|5.13|) again, we deduce 

|T 2 |<^M(/,5)KHo,n. (5.19) 
Taking ([5TT8]) and (pTlUj) into (|5TT7|) . the desired result (f37TU]> is obtained. □ 

Finally we give the error estimate of the trace flux . 
Theorem 5.3. Let A^ be the solution of (|5.1f) - (|5.4|) . Then 

\\u-\ h \\ , d T h <^M(f,~g). (5.20) 



Proof. For e G &Th\dVL, 



Ah = {{t4}} + ^[QJ- 



Note that q ■ n is continuous across inner edges(faces). By Lemmas I4.3H4.4I we have 

\\ u ~ ^h\\o,ar h \an < Ik - ^/»llo,a7^\an + tHIII? ~ QJIIo,97i\<3^ 

< ||u - n h u|| 0i9Th + ||n h u - ^||o,9Th + t -1 (II<? - n-hq\\o,8T h + \\ u hq - QJksrJ 

< \\u - IL h u\\ 0t dT h + ph~*(\\u - U h u\\ 0> n + \\u - U h \\ ,n) + - H-hq\\o,dT h 

+ r- 1 p7»"'(l|9-n hq || 0l n + ||g-Qhl|o,n) < —M(f,g), 

p 
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where we utilize Theorems I5.HI5T21 in the last inequality. 
For eedT h n dQ, 

t(u -U h ) (q-Q h )-n g- P M g 

u ~ X h = — j— + 7- h — j— ■ 

1 + r 1 + r 1 + r 

Since 

ll« - X h\\l,an = ll u - n Hlo,an -\\Xh~ U h u\\ 2 dn - 2Re(u - X h , X h - U h u) dn , 
the definition of Pm implies that 

2r 

\\ u ~ ^h\\o,aQ < \\ u ~ ~Rhu\\l.dn ~ W^h - H-hu\\l dn + — — \\u - U h \\ ,du\\Xh - n fe n|| ,an 
2 

+ t~ ; — — Q/illo,9nll^/i _ n^wLan 
1 + r 

< [|u - n ft «||g >an + ||u - J7h||o,an + T ~ 2 \\l ~ Qh\\o,m- 
Similar to the proof for the case of inner edges, we can also show that 

||«-A A ||o,ai < —M(f,g). 
p 

Hence (|5.20p is derived. □ 

6 The error estimates 



In this section, we shall derive error estimates for the scheme (|2.5p - (|2.8|) . This will be done by 
making use of the stability estimates derived in Theorem l4.1l and the error estimates of the auxiliary 
problem established in the previous section. 

Theorem 6.1. Let q h , Uh and Uh be the solution of (|2.5p - (|2.8|) . We have 

. /nfi 2 K 2 h 2 k 5 /i 4 \ . ~ , , 
«-«fc|o,n< + — 5- + — S- 



V p z p z p* , 

. fnh n 3 h 2 k 6 /i 4 \ , ~ . , 

4Q - Qhhfl ^{— + — + — W(f,g), 6.2 
\ p p z p^ J 

^ /tzhi ,Kh,S. ,nh,5.K 3 h 2 K?hh K 5 h^\ ,~ . . 

u-u h , dTh < + — 2 + — ) 2 — + H —)M(f,g)- 6.3 

\ p p ' p p z p p 6 J 



Moreover, if < 1, then 



||«-«h||o,n< (^ + ^)M(/,2), (6.4) 

4Q-Q h \\o,n<(- + ^PjM(f,g), (6.5) 
\ p p L / 

\\u-u h \\o,dT h < + + s-)M(f,g). (6.6) 

V p p p 2 ' 
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Proof. Denote := q h — Q h , e\ := Uh — Uh, and := Uh — A/,, according to the formulation 
H223I-II23D and dSHD-<EI3D, we have (e£,e£,ef) G F£ x x and they satisfy 



(ine q h ,r) Th - (e£,divr) Th + (e u h: T^n) dTh = 0, 

^)r h - (e*, Vw) Th + (q h ■ n,w)er h = 2\K{U h - u,w) Th , 
(-^•n + ef,/J}an = 0, 
{th-n,~P)dT h \dn = 0, 
for all r G V?, G VF?, and /i G Mf. The numerical flux is given by 

i q h = 4 + r(e u h -ei)n on dT h . (6.7) 
The stability estimates in Theorem 14. 1 1 imply that 

(Kj JlP \ 

(Ki*^ hp 1 \ 
1 + ^p-) K h-U h \\ 0t n, 

\\e%,ar h < ((y) 1 +ph~*) (l + ^-)«||« - U h \\o,a. 
Using the triangle inequality and Theorem 15.21 the proof is finished. □ 
Next we demonstrate the improved convergence results for the coarse meshes under the condition 

3 1,2 

> 1. First we give the stability estimate for the following elliptic HDG scheme. 

Lemma 6.1. Let (Q h , U h , U h ) G V\ x x M p be the solution of the following elliptic HDG 
scheme 



{wQ h ,r) Th - (U h ,divr) Th + (U h ,r^n)gr h = 0, 
-{\tdA h ,w) Th - {Q h ,Vw) Th + (Q h ■ n,w) dTh = (f,w)r h , 

{-Qh • n + U h ,~p)dn = 0, 
(Qh • n,Jl) dTh \ dn = 0, 

for all r G V p h , w G W^, and fj, G M^, where the numerical flux Q/i is given by 

Qh = Qh + r(Uh - U h )rt on d%. 

Then there holds 

||Qft||o,n + ||Z4|| 0l n<«- 1 ||/||o,n. (6- 
Proof. Similar to the proof of Lemma 14.51 we can deduce 



{\nU h Mh)T h + (usQh, Qh)r h + {r(U h -U h ), {U h -U h )) dTh + {U h ,Uh)m = (fM%- 

Hence 

«||W ft ||^<||/||o,n||^|| 0j n, 
«l|GAllo,n < ||/IM|Z4||o,n, 

which means 

K(||Z4|| 0l n + ||Qh||o I n) < ||/||o,n. 



□ 
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Theorem 6.2. Let q h , Uh and be the solution of (|2.5p - (|2.8p . Under the mesh condition > 1, 
we have 

||u-^|| 0) n<^M(/,5), (6.9) 

„ (Kh K 3 h 2 \ , ~ „. , 

K 9 " fto,!! ~ 1 —)M(f, g), 6.10 

\ p p z ) 

Il« - fi/»||o,ar h < + —r)M{f,g). (6.11) 

Proof. In this mesh condition, the stability estimates in Theorem 14.11 and Lemma 14.51 indicate the 
following inequality 

3 

hh-u h \\ ,ffr h <—M(f,g). (6.12) 

p2 

The consistence of the HDG scheme implies that 



- Qh),r)T h ~ (u - u h ,divr) Th + (u- u h ,r^n) dTh = 0, (6.13) 

(ik(u - u h ),w) Th + (div (q - q h ),w) Th - r(u h - u h ,w) dTh = 0, (6.14) 

(-(<? - Qh) • n + ( u ~ Uh),~P)an = 0, (6.15) 

((Q-Qh)- n ,~P)dT h \dSl = 0, (6.16) 

for all r E V?, if G W?, and /u G Mf. We introduce the dual problem which replaces the right 
hand side of (|4.2ip - (|4.23p by u — Uh as follows: 

-i/e* + Vtf = infl, (6.17) 
div * — in^ = u — Uh in $7, (6.18) 
* • n = * ondn. (6.19) 

Similar to Lemma |4.6| we have the following regularity estimate: 

K~ 2 \\^h,n + « _1 [|¥[|i )n + « -1 ||*||i jn < ||u - «fc||o,n. (6.20) 
We denote by Pm the I? projection onto M£, 

(P M W, V)dT h = tAdTh for all G M£. 
From (|6.13p and (|6.18p . we can easily get 



(u -u h ,u- u h ) Th = (u - u h , div * - iK^f)r h = {u- u h , div $) Th + ik(u - u h , 



= (u - u h , div (* - n h *)) Th + (i«(g - q h ),n h &) Th + (u-u h ,U h ^ ■ n) dTh +in(u - u h ,^) Th . 
According to (|6.17p and Green formulation, there holds 

(HQ ~ Qh),*)T h = -(Q - Qh,^)% = (div (q - q h ),^) Th ~ ((Q ~ Qh) ' n^)dT h - 
Hence, by (|6.14p we obtain 

(u -u h ,u- u h )r h = Ai + A 2 , (6.21) 
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where 

Ax = (u-u h ,div($ -Il h $>)) Th - (iK(q - q h ), & -Il h &) Th - {{q - q h ) ■ n,V - P M ^)dT h 
+ (i/s(u - u h ) + div (q - q h ), ,V -Il h V) Th - t(u - u h ,U h ^ - P M ^)dT h , 

and 

M = r(u h - u h , n h y)dT h ~ ((<? ~ Qh) • n - ^)&T h +{u- u h , Il h $ ■ n) aTh 
+ ((q - q h ) ■ n, * - P M ^)dT h + r(u - u h , U h ^ - P M ^)dT h - 
Utilizing (|6"7T5|) . ([6361) and (l6TT9|) . A 2 becomes 

M = -((<? - Qh) ■ n > P M^)dT h + T (uh ~ Uh, PM^)dT h + T ( u h ~ U h ,U h ^ - P M ^)dT h 
+ {u- u h , (Il h & - • n) dTh + {u- u h ,W- n) m + t{u - u h ,Yl h ty - P M ^) dTh 
= -{u- u h ,P M ^) d n + t{u - u h ,n h ^ - P M V}dT h 
+ (u- u h , (II h & - *) • n)ar h + (u-u h ,W- n) d n 

= (u - %, * - Pntyan + r(u - u h ,U h ^ - P M *)dT H + ( u ~ «f» ( u h® ~ *) • n) dTh . 

Taking the complex conjugation of (|6.21|) and making use of Green formulation and the definition 
of projection operator, we obtain 

h - «fc||o,n = - U h ^,q^q^) Th + (* - n h *,div (q - q h )) Th - (* - Pm*, (q ~ Qh) • n) dTh 

- - n h v,u-u h ) Th - (* - n / ,*,v(u-u/ l )) Th + <(* - n h *) • n,u - u h ) dTh 

- r{n h ^ - p m $,u- u h ) dTh + (* - p A /^,n-u h ) 9c 

- ((* - n^*) • n, u - u h ) dTh + T{U h ^> - P M V, u - u h ) dTh 

= - n h *,g- n^)^ + (* - n^div q) Th - - p M #, (g - n^g) • n)ar h 
-i/t(*-n^,u-n ft «) Th - (*-n A *,v(u-7r£u))r h - ((*-![**) ■n,« fc -tt A ) arh 
+ r(n ft tf - Pm*,^ - u/Oar* + (* - Pm*, « - n fc «)an- 

Note that 

||n h ^ - P W ¥||o,Sr h < ||* - n h *||o,9T h + II* - PM%\o,dT h < II* " W\\o,dT h - 
Using Lemma |4~4"1 the regularity estimates (|2.4p and (|6.20p . the following inequality is derived 

K 2 h 2 

\\ u ~ u h\\o,n ^ — 2~\\ u ~ Uh\\o,uM(f,g), 

which means 

K 2 h 2 

\\u - u h \\ ,Q < — =-M(/, 5). 

Since (e^,e^,e^) G V p h x x satisfy 

( iKe Mr h - (eh,divr) Th + {ef l ,r^n) dTh = 0, 
"Ki^Ti - i e t Vw )T h + (Qh ■ n MdT h = 2i/s(u - u h ,w) Th , 

(-£ Q h -n + 4^)dn = 0, 
<efc-».M)s7fc\an =0, 
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for all r E V^, w E and /i E Mf, it follows from the stability estimate in Lemma 16. II that 
K \\q - Qh\\o,n < it\\q - Q h \\o,n + «[|e)[[|o,fi < ( 1 —)M(f,g). 



p jr 



The triangle inequality and (|6.12|) imply that 



W - Uh\\o,dT h < 11^ - H-hu\\o,dT h + \\ u h - n^nllo,^ + \\ u h - Uh\\o,dT h 

-1 h 5 

< ph~{\\u - n fc u||o,n + \W - u h \\ ,n) + K(-) 2 M{f,g) 



-2, 3 - 3 



< 



K z fi2 nhi ^ , ~ „. 



P p2 

The proof is completed. □ 

7 Numerical results 

In this section, we present a detailed documentation of numerical results of the HDG method for 
the following 2-d Helmholtz problem: 

-Au-A = /:=^ infi, (7.1) 

dn 



^+\KU = g on T R := dQ. (7.2) 



Here O is unit square [—0.5,0.5] x [—0.5,0.5], and g is chosen such that the exact solution is given 
by 

cos«;r cosK + isinK , 

U = . ; . r , 7 , ,, Jo(Kr) (7.3) 

k k(J q (k) +iJi(«)) 

in polar coordinates, where Jjy(z) are Bessel functions of the first kind. 

In the numerical results of [T7j, the optimal convergence of the HDG method is observed when 
the parameter r is chosen as 0(1). In this work, when u S H 2 (£l) we let r = which is also used 
in the following experiment. The HDG method is implemented for piecewise linear (HDG-P1), 
piecewise quadratic (HDG-P2) and piecewise cubic (HDG-P3) finite element spaces. 

For the fixed wave number k, we first show the dependence of the convergence of ||« — «h||o,fi, 
\\q~ Qh\\o,n an d \\u — Uh\\o,8T h on polynomial order p and mesh size h. On one hand, the left graphs 
of Figure [2] display the above three kinds of errors for k = 100 by HDG-P1, HDG-P2 and HDG-P3 
approximations. We find that the pollution errors always appear on the coarse meshes, but the 
errors of \\u — Uft||o,n almost converges in 0(nh 2 /p 2 ) on the fine meshes, and \\u — Uh\\o,dT h nearly 

3 

converges in 0(nh2 /p) on the fine meshes. The results support the theoretical analysis. We note 
that the error of ||qf — Q^||o,Q a ls° almost converges in 0(k/i 2 /p 2 ) on the fine meshes, which is a little 
better than our theoretical prediction. On the other hand, for the case of k = 300, the right graphs 
of Figure [2] show that the errors of \\u — Uh\\o,dT h i \\ u ~ u h\\o,n and \\q — gJ|o,n always decrease for 
high order polynomial approximations. 

Figure [3] displays the surface plots of the imaginary parts of the HDG-P1, HDG-P2, HDG-P3 
solutions of Uh and the exact solution for k = 100 with mesh size h ~ 0.022. It is shown that the 
HDG-P2 and HDG-P3 solutions have correct shapes and amplitudes as the exact solution, while 
the HDG-P1 solution has a correct shape but its amplitude is not very accurate near the center of 
the domain. 
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■ llq - q l. .„ 

- L2 -error of u on the trace 

■ 0(K h 2 /p 2 ) 

■ 0(K h M /p) 



10 



HDG-P3(K=100) 








Figure 2: Errors of ||u — Uft||o,n, \\q — <Z/Jo,n and \\u — Uh\\o,dT h f° r K = 100 by HDG-P1, HDG-P2 and 
HDG-P3 approximations (left, top bottom); Errors of \\u — Uh\\o,dT h > — u h\\o,n and ||q — <7/J|o,n for n = 300 
by different polynomial approximations (right, top bottom). 
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HDG-PI Solution of u h HDG-P2 Solution of 




Figure 3: Surface plot of the imaginary parts of the HDG-PI, HDG-P2, HDG-P3 solutions of Uh and the 
exact solution for k = 100 with mesh size h s» 0.022. 




50 100 150 200 5 10 15 20 25 30 35 40 45 50 55 



Figure 4: Errors of \\u - u h \\o,dT h ■ Left: ^ = 1.1. Right: = 1.1. 
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Figure 6: The error of \\u - u h \\ Q . aTh for re = 50, 100, 200, 300 by HDG-P1, HDG-P2 and HDG-P3. 
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The left graph of Figure 0] displays the error \\u — Uh\\o,dT h f° r fixed ^ = 1.1. It shows that the 
error \\u — Ufc||o,97^ can not be controlled by ^ and increases with k, which indicates that there 
is pollution error in the total error. The right graph of Figure H] displays the same error with the 

3 u2 

mesh size satisfying = 1.1 for different k, h and p. We observe that under this mesh condition, 
the error \\u — Uh\\o,dT h does not increase with k. The top two graphs of Figure [5] shows the similar 
property for the error n\\q — Q/J|o,Q. From the bottom graph of Figure [5] we can also find that the 
error \\u — Uh\\o,Q does not increase with k only under the mesh condition ^ = 1.1. 

Next we verify the convergence properties of the HDG method for different wave numbers by 
piecewise PI, P2 and P3 approximations respectively. In FigureEl the error \\u — Uh\\od%, of HDG- 
Pl, HDG-P2 and HDG-P3 solutions for different wave numbers always oscillates on the coarse 
meshes, and then decays on fine meshes. For HDG-P1 solution, \\u — Uh\\o,dT h grows with k along 
the line nh = 0.8 for k < 300. By contrast, for k < 300, this error does not increase significantly 
along the line nh = 0.8 by HDG-P2, and decreases along the line nh = 0.8 by HDG-P3. This also 
means that the pollution error can be reduced by high order polynomial approximations. Figure 
[7] shows the convergence property of \\u — li^ ||o,Q f° r different wave numbers. For k < 300, along 
the line nh = 0.8, the error \\u — Uh\\o,n of HDG-P1 and HDG-P2 solutions stays stable, and for 
HDG-P3 solution, this error decreases as \\u — Uh\\o t dT h - Similar phenomenon can also be observed 
for the error [|<f — <f/J|on by different polynomial approximations. 

HDG-P1 HDG-P2 




1/h 1/h 



HDG-P3 

0" 1 f ■ ■ — . . — — 




10 ' ' — ' ' — ' ' — 

10° 10 1 10 2 10 3 

1/h 

Figure 7: The error of \\u - u h \\ ,n for k = 50, 100, 200, 300 by HDG-P1, HDG-P2 and HDG-P3. 

For more detailed comparison between HDG methods with different polynomial order approx- 
imations. We consider the problem with wave number k = 200. The traces of imaginary part of 
the HDG solution with piecewise PI, P2 and P3 approximations in the rcz-plane with mesh 
sizes h ~ 0.022, 0.0055, and the trace of imaginary part of the exact solution, are both plotted in 
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Figure [HI On the coarse mesh with h ~ 0.022, the shapes of HDG-P2 and HDG-P3 solutions are 
roughly the same as the exact solution, while the shape of HDG-P1 solution does not match the 
exact solution well. But on the fine mesh with h w 0.0055, the shapes of HDG solutions match well 
and even better for high order polynomial approximations. Then we can observe that although the 
phase error appears in the case of coarse meshes and low order polynomial approximation, it can 
be reduced in the fine meshes or by high order polynomial approximations. 
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Figure 8: The traces of imaginary part of the HDG solution u h by HDG-P1, HDG-P2 and HDG-P3 (top 
downbottom) with mesh sizes h w 0.022,0.0055 (left, right). The trace of imaginary part of the exact 
solution is plotted in the green lines. 
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